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SUMMARY 

Fouling (and corrosion) of gas turbine (GT) blading often results from 
the deposition of materials derived from Inorganic Impurities In the fuel 
and/or Ingested air. When the depositing material Is In the form of a sub- 
micron "aerosol" (dust or mist) It has been found that the rates of deposition 
on cooled surfaces can be augmented by some 100-1000-fold via the mechanism of 
thermophoresis (particle migration down a temperature gradient). For this 
reason. In earlier papers we reported the development of rational, yet simple, 
engineering correlations of thermophoretlcal ly-augmented particle transport 
across both laminar boundary layers (LBLs) and turbulent boundary layers 
(TBLs) . While developed based on theoretical considerations, and numerical 
computations of self-similar LBLs and law-of-the-wall (Couette flow-like) 

TBLs, these mass transfer coefficient (St m -) correlations, when applied 
locally, may also prove useful In making engineering mass transfer predictions 
for more complex geometries, Including GT-blades. Pending additional con- 
trolled experiments, Insight Into the local applicability of these cor- 
relations can be gained by selected comparisons with numerical predictions for 
developing BLs. This paper reports on our (a) modification of the code STAN5 
to properly Include thermophoretlc mass transport, and (b) examination of 
selected test cases of developing BLs which Include variable properties, 
viscous dissipation, transition to turbulence and transpiration cooling. 

Under conditions representative of current and projected GT operation, local 
application of our St m /St m 0 - correlations evidently provides accurate 
and economical engineering design predictions, especially for suspended 
particles characterized by Schmidt numbers outside of the heavy vapor range 
(say, Sc>10). 
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INTRODUCTION 


Accurate predictions of mass transport rates In nonlsothermal forced 
convection systems are now essential to the gas turbine (GT) Industry, among 
other technologies. Recently, a dramatic enhancement In small particle d 1 f- 
fuslonal transport rates due to thermophoresis (particle drift down a temper- 
ature gradient) has been predicted for both Internally cooledand trans- 
piration-cooled surfaces across both self-similar laminar boundary layers 
( LBLs) (refs. 1 and 2) and law-of-the-wall turbulen t (ref. 3.) boundary layers 
(TBLs), Including viscous dissipation (ref. 4). The magnitude and temperature 
dependence of this thermophoretlc augmentation has been experimentally 
confirmed for LBL flow of combustion products containing submicron 
MgO-partlcles (ref. 5). 

To facilitate engineering predictions, simple rational correlations have 
recently been developed (refs. 6 and 7). Until now their accuracy and 
applicability limits have been demonstrated using “exact" numerical 
calculations for self-similar BLs (refs. 1 and 4). One of our purposes here 
Is to examine their behavior and accuracy In more complex BL-flows, under 
conditions representative of current and future GT-technology . 

The availability of high-speed, large capacity computers has made It 
possible to examine the reliability of self-similar heat transfer correlations 
when applied locally to developing BL situations. Indeed, computer codes 
for both multidimensional laminar and turbulent flows have been developed by 
exploiting finite-difference techniques (refs. 8 and 9). Recently, with 
Increased Interest In GT blade fouling and/or corrosion problems, greater 
attention Is being focused on small particle mass transfer across two- 
dimensional BL 1 s (ref. 10). To make such predictions we have adapted the 
two-dimensional BL code STAN5 (ref. 11), which has previously been widely 
used for gas-side convective heat transfer predictions. This necessitated 
modifying the program to now Include thermal (Soret) diffusion 
( "thermophoresls" for small particles) In the suspended particle mass con- 
servation (particle-phase "continuity") equation. Toward this end, STAN5 
predictions over the suction surface of a typical stator blade are here 
compared with the predictions of our earlier correlations (refs. 6 and 7) when 
applied locally to such developing BL situations. Of course, experimental 
data, both for blade cascades and operational turbines, will be the ultimate 
arbiter of whether our correlations and/or numerical BL codes (e.g., STAN5) 
presently and properly Include all of the phenomena which dictate mass 
transport rates In the challenging GT environment. 

NOMENCLATURE 

B^ real blowing parameter, eq. (5) 

-B-j. thermophoretlc suction parameter, eq. (6) 

c mass fraction of particles In the local gas mixture 

D Brownian diffusion coefficient of particles 

Da effective Damkohler number, eq. (7) 

Le Lewis number (ratio of particle Brownian dlffuslvlty to host gas 

heat dlffuslvlty) 
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m" mass flux at station e or w 

Ma Mach number 

p pressure 

r radius (for axlsymmetrlc BL flow) 

r nose radius of the turbine blade 

n 

S source; eq. (3) 

Sc Schmidt number (ratio of host gas momentum dlffuslvlty (kinematic 

viscosity) to particle Brownian (mass) dlffuslvlty) 

St. Stanton number for heat transfer 

h 

St Stanton number for mass transfer, Including thermophoresis 

m 

St Stanton number for mass transfer without thermophoresis and/or 

m,o 

transpiration cooling 
T absolute temperature 

T temperature at the outer edge of Brownian diffusion boundary layer 

m 

u fluid velocity In x-dlrectlon (parallel to wall) 

v fluid velocity In y-dlrectlon (normal to wall) 

x distance along the surface (measured from the forward stagnation 

point) 

y distance normal to surface 

thermal diffusion factor of the particle/host gas combination 
p dynamic viscosity of host gas 

p density of host gas 

y stream function coordinate, pur = 3\|//ay, pvr = 

« nondlmenslonal stream function coordinate, (^'-'J' w )/('l» e -'l» w ) 

Subscripts: 

e pertaining to the outer edge of the boundary layer 

eff effective 

o pertaining to reservoir conditions (transpiration coolant) 

w pertaining to the surface (wall) 

oo pertaining to gas mixture at upstream "Infinity" 

Miscellaneous : 

BL boundary layer (L, laminar; T, turbulent) 

GT gas turbine 

PDE partial differential equation 

RHS right hand side 
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ADDITION OF THERMOPHORESIS TO STAN5 BL CODE 


Here we focus only on the suspended particle mass conservation POE, 
which must Include the thermophoretlc flux term. Upon time-averaging, the 
simultaneous fluctuation of the thermophoretlc speed and particle concen- 
tration produces a new "correlation" term. But this latter term (which may be 
called "eddy thermophoresis") Is neglected based on arguments given In 
references 12 and 13. In the same notation as reference 11, the time-averaged 
particle mass conservation equation for the BL can then be written: 


3c_ 

ax 


ac 

ay 



( 1 ) 


where the second term on the RHS of equation (1) Is the added thermophoretlc 
flux. Upon Introducing the nondlmenslonal stream function as the transverse 
coordinate (cf. ref. 11) equation (1) becomes: 
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In the numerical solution of equation (2) the thermophoretlc flux term on 
the RHS will be treated as a "source" term. Although the source Itself Is 
explicitly c-dependent, c-values obtained from the previous (upstream) 
Integration step are used In Its evaluation. This simple treatment of the 
thermophoretlc flux term Is compatible with numerical approximations already 
Incorporated In STAN5. For flow near the wall, where STAN5 switches to 
simplified Couette-flow equations (with no streamwlse variation of the 
dependent variable, and constant r), the thermophoretlc flux term Is Included 
In the particle mass conservation equation, again as the "source" defined by: 


S 



(3) 


Based on the procedures outlined In reference 14, new subroutines have 
been added to STAN5 In order to calculate the suspended particle thermal 
diffusion factor, oy, and the Brownian diffusion coefficient, D (as well 
as the Schmidt number. Sc), Including their respective temperature 
dependencies . 
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The modification In particle disposition rate due to thermophoresis 
depends both on particle size (hence Schmidt number) and the ratio of 
malnstream-to-surface-temperature. Therefore, under typical turbine operating 
conditions, In order to see an appreciable thermophoretlc effect one needs to 
consider particles heavier than ordinary vapor molecules. However, as Is well 
known, when the particle size Increases (enough to cause Sc » 1), the mass 
transfer (Brownian diffusion) BL Is completely Imbedded within the momentum 
(or heat) transfer BL, the thickness ratio being proportional to Sc -1 ^ 3 
(or Le + '^ 3 ). In fact. It can be readily shown (ref. 1) that In the 
presence of "strong" thermophoresis ( 1 . e . , aj [(T e /T w ) - 1] » 1) the 
Brownian diffusion BL Is even much thinner than In the absence of 
thermophoresis, with the thickness ratio scaling becoming proportional to 
Sc -1 (or Le + '). These observations Imply that a finite-difference 
scheme with a much finer mesh size Is required for an accurate solution of the 
small particle concentration profile and the associated mass transfer Stanton 
number. Moreover, to preserve accuracy for a two-dimensional Integration 
process, a thinner BL requires a smaller step-size In the x-dlrectlon — 

1 .e. . x-dlrectlon step-size Is restricted by the smallest of the momentum, 
heat and mass transfer BL thicknesses. Consequently, a moderate Increase In 
the particle size (via Sc) Increases the computation time dramatically. In 
order to achieve a reasonable balance between (a) sufficient thermophoretlc 
deposition rate enhancement (to test the accuracy of our correlations) and (b) 
required computation time, we have therefore adopted a Schmidt number of about 
26 for most of our trials (except for one case with Sc = 100, see "Results 
and Discussion"). Accordingly, the dimensioning of STAN5 has been Increased 
up to 1000 mesh points In the y-dlrectlon (from the originally programmed 
maximum of 60 mesh points) 4 . 


RESULTS AND DISCUSSION 

For all results presented here a mainstream velocity distribution 
obtained from reference 15 for a turbine stator blade Is used as Input to 
STAN5. Except for figure 6, the nose radius, R n , of the blade used for the 
illustrative calculations was taken to be 0.152 cm. Integrations were carried 
out along the suction surface up to about 10 nose radii from the blade stag- 
nation point. In these cases the predicted BL flow remained laminar. To 
examine the consequence of transition to turbulence (by x/R n «5), for 
figure 6 the blade was magnified to R n = 0.366 cm and the stagnation 
pressure was Increased from 4 atm. to 20 atm. 

In all figures, results at the stagnation point (x = 0) are obtained from 
the self-similar solutions of reference 1 (because STAN5 is singular at x = 

0) . STAN5 calculations are started at x = 10 _1 R n . Downstream results 
are very sensitive to the Initial starting profile specified (see also 


4 For the fastest commercially available computer, CRAY-1, using 800 
mesh points In the y-dlrectlon and taking the step-size along the x-dlrectlon 
to be 10"^ times the local momentum BL thickness, Integration of LBL 
equations by STAN5 along the suction surface of a typical turbine blade to 
about 10 nose radii downstream from the stagnation point required a CPU time 
of about 8 minutes with thermophoresis and 4 minutes without thermophoresis. 
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ref. 16). To be on the safe side, only STAN5-values beyond one nose radius 
from the stagnation point were considered reliable (Initial profile In- 
sensitive). Moreover, by moving one nose radius away from the stagnation 
point we avoid any anomalltles that may result due to (a) lack of accuracy In 
the mainstream velocity near the stagnation point (ref. 15), and (b) an un- 
realistic curve-fit generated by the cubic splines that STAN5 employs (ref. 
17). The automatic starting profile generating subroutine of STAN5 (ref. 16) 
was used for the Initial velocity and temperature profiles, and the velocity 
profile was scaled (via Sc -1 / 3 ) to obtain the Initial concentration profile. 

The correlation values shown In all figures (with negligible viscous 
dissipation) are calculated from (refs. 6 and 7): 


where: 
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The Schmidt numbers reported In the figures correspond to particles 
(assumed spherical) made up of Na 2 S 04 "monomer units" In air 5 . The 
ratio ajLe (of the thermophoretlc dlffuslvlty of the particles to the heat 
dlffuslvlty of air), required for our correlations (ref. 6) was estimated to 
be 0.4 and temperature-independent for all reported cases here. For 
simplicity In the correlation scheme, throughout this paper the ratio: 
Sth/St m>0 , of the Stanton number for heat transfer to the Stanton number 
for mass transfer without thermophoresis, was approximated to be equal to the 
prevailing value of Le w -2/f3 for laminar BL flows. However, for 
transitional and fully turbulent flows, local St^ and St m 0 values 
(as calculated by STAN5) are used for the correlation calculations. 

To "calibrate" STAN5, we first ran the case of a flat plate (constant 
mainstream velocity), with the results shown In figure 1. The abscissa Is the 
distance along the flat plate made nondlmenslonal by the stator nose radius 
R n = 0.152 cm to allow comparisons with the abovementloned turbine blade 
example. The ordinate Is the ratio of mass transfer Stanton numbers 


5 Sc CO - 26 corresponds to 10^ and Sc = 100 corresponds to 8x10^ 
Na 2 S 04 monomer units In a cluster. 


6 


(obtained with thermophoresis and without thermophoresis, respectively). 

Under the prevailing conditions there Is more than a two-fold enhancement In 
the deposition rate due to thermophoresis. For this "self-similar" case the 
result obtained from the exact numerical calculations of reference 1 Is shown 
by the dash-dot line (x-lndependent) . The correction-predicted value Is the 
dashed line. The only difference between the results of reference 1 and STAN5 
should be that due to the somewhat different treatment of variable host gas 
(air) thermodynamic and particle transport properties for this "cold" wall 
(Tw/Tc = 0.6) self-similar case 6 . However, In Interpreting subsequent 
comparisons, the reader should recall that STAN5 does not quite give a 
constant (x-lndependent) enhancement. 

Figure 2 displays the results of calculations for the suction surface of 
a stator blade, using the same environmental conditions as In figure 1. For 
the prevailing conditions the Mach number Is calculated to be about Ma e = 

0.72 at x/R n = 10 (less for smaller x/R n ). To display the effect of 
the Mach number STAN5 was also run Including viscous dissipation. As expected 
(ref. 4), at these Mach numbers viscous dissipation Is seen to have a small 
effect on the thermophoretlc enhancement of deposition rates. Comparison of 
our correlations with STAN5 both with and without viscous dissipation shows 
that the local difference Is always less than 10 percent. The shape of the 
curves (l^e. , the fact that the effect of thermophoresis decreases as x/R n 
Increases) can be explained by the mainstream static temperature cooling. 
Although T w /T,„ Is fixed for our calculations, the difference between 
the mainstream static temperature, T e (which gets cooler as the flow 
accelerates under adiabatic conditions), and the wall temperature, T w 
(constant), gets smaller, which reduces the effect of thermophoresis. Of 
course, our correlations are also applied locally, with the relevant local 
temperature being T e rather than T,,,. 

Figure 3 Is constructed for the same conditions as figure 2 except that 
the wall temperature Is higher (Ty/T* = 0.8), corresponding to a smaller 
enhancement In deposition rate due to thermophoresis. The agreement between 
the predictions of STAN5 and our correlations Is excellent. 

In figure 4 an Increased particle size Is considered (Sc,*, = 100), 
keeping all remaining parameters the same as In figure 3. The corresponding 
effect of thermophoresis Is now even greater than In the colder wall case 
(Ty/Too = 0.6) with Sc® = 26 (shown earlier In fig. 2). Yet, the 
agreement between our correlations and STAN5-pred1ct1ons Is even better, with 
local differences being less than 4 percent for all x/R n -values. This may 
be related to our observation (ref. 6) that for self-similar cases these 
correlations Improve as T w /Ta, -♦ 1 , regardless of particle size (Sc). 

To consider even larger departures from self-similarity a variable 
(x-dependent) wall temperature was considered, with the results shown In 
figure 5. T w Is varied In such a way that the wall Is hottest at the 
stagnation point (Ty/T*, = 0.8) and gets colder as x/R n Increases, 
becoming Ty/T*, = 0.5 at x/R n = 10. Not surprisingly, the 
thermophoretlc effect Increases as the wall gets cooler (although the 


Reference 1 assumes that the host gas thermodynamic and particle 
transport properties have simple power-law temperature dependencies. 
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mainstream static temperature drops as the flow accelerates, the reduction In 
the assumed wall temperature Is greater, resulting In a steeper temperature 
gradient across the BL) . While departures between the correlation and 
STAN5-pred1ct1ons Increase as the wall-to-malnstream-temperature ratio departs 
from unity (ref. 16), remarkably enough, the effects due to "nonself- 
slmllarlty" act In such a way that the agreement between STAN5 and the 
correlations Improve as T w drops. 

In figure 6 the nose radius of the stator blade and the pressure have 
been Increased to 0.366 cm and 20 atm, respectively. The corresponding BL 
flow becomes transitional at about x/R n = 2.9 and "fully" turbulent at 
about x/R n = 5.0. A lower wall temperature (T w /Too = 0.5) Is considered 
to provide a stringent test of the present correlations, based (ref. 6) on 
self-similar Integrations for T w /Too > 0.5. We observe that the 
agreement between STAN5 and our correlations Improves for turbulent BLs. In 
this connection It should be recalled that our correlations describe the 
results of law-of-the-wall Couette-flow TBL Integrations with even better 
accuracy than In otherwise corresponding LBL situations. When compared 
locally with STAN5 results, even for T w /T,a = 0.5, agreement Is within 
10 percent for the laminar portion, and within 4 percent for the fully 
turbulent portion of the flow. The agreement Is still better In the 
transitional region. Note that the effect of thermophoresis on mass transfer 
(over and above the already efficient transfer mechanism by turbulent eddies) 
Is smaller than In the corresponding laminar BL case. 

Figure 7 displays the results of our calculations on the effect of 
transpiration cooling and/or thermophoresis on the deposition rate of small 
particles (Sc,,, = 26). For this purpose STAN5 was run with the wall 
boundary condition specified as a flux (Instead of a level), where the 
reservoir temperature of the coolant air was taken to be 600K (Tq/To, = 

0.4). The correlation curve without thermophoresis Is representative of how 
the blowing rate was varied along the blade surface. The blowing rate was set 
higher near the stagnation point, was reduced along the suction surface, and 
near the trailing edge (7 < x/R n < 10) was chosen such that the real blowing 
parameter, B m (defined by eq. (5)), was constant. 

If attention Is first focused on the STAN5 and correlation curves without 
thermophoresis, observe that STAN5-pred1ct1ons "lag" In responding to local 
blowing rate variations when compared with the correlation (which Is based on 
an "Instant" assumed adjustment to the local blowing rate). Specifically, 
for 7 < x/R n < 10, although the correlation predicts a constant reduction 
In the deposition rate, STAN5-values steadily Increase, ultimately approaching 
the correlation prediction, but remaining below them due to "recollection" of 
the "history" of higher upstream blowing rates. For a two-dimensional 
developing BL each x-lntegratlon station provides the starting profile for 
the next station In the forward "marching" process. Therefore, the flow 
cannot Instantaneously adapt Itself to the “new" (continuously varying) 
boundary conditions. A better example Is provided If the cross-over points of 
the curves with and without thermophoresis are compared for the correlation 
(x/R n «7.4) and STAN5 (x/R n «8.2). Because of mainstream static temper- 
ature cooling and the small blowing rates used In this example, the wall 
temperature becomes hotter than the mainstream static temperature at about 
x/R n «7.4. Although thermophoresis enhances small particles deposition 
rates toward colder walls ("suction" effect, ref. 6), for hot walls the 
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deposition rate Is reduced due to what we have called "thermophoretlc 
blowing." The correlation curve with thermophoresis, of course, exhibits an 
Immediate reaction to the hotter wall, and crosses the correlation curve 
without thermophoresis at x/R n ~7.4 (T w /T e = 1). However, STAN5-w1th 
thermophoresis responds to the change to a hotter wall further downstream at 
x»8.2 and predicts lower deposition rates than STAN5 - without thermophoresis 
from then on. We estimate the "lag" In the reaction of STAN5 to be about 9 to 
10 local momentum BL thickness for the conditions Investigated here. For a 
blowing rate distribution with more abrupt changes we can therefore expect 
even larger local disparities between STAN5-pred1ct1ons and such cor- 
relations. However, In such cases (with sharp boundary condition variations) 
It should be recalled that the reliability of STAN5 results themselves becomes 
questionable 7 . Considering the cost and accuracy features of both methods, 
the comparisons shown In figure 7 suggest that local applications of our cor- 
relations will be satisfactory for most engineering design predictions on 
transpiration-cooled objects. 


CONCLUDING REMARKS 

To examine the behavior and accuracy of our recently proposed 
St m /St m ^-correlations of particle mass transport (ref. 6) when applied 
locally’to developing GT-boundary layer situations, we have embarked on a 
program of selected numerical computations and controlled laboratory 
experiments (ref. 5). Progress In the former category Is reported here, and 
In view of the widespread familiarity and use of the two-dimensional BL code 
STAN5 (ref. 11 and 9) for GT heat transfer predictions, we have described 
how STAN5 has been adapted for mass transfer (deposition) rate predictions for 
cooled combustion turbine blades. Most Importantly, the program has been 
modified to Include thermal (Soret) diffusion ("thermophoresis" for small 
particles) In the suspended particle mass conservation equation. The 
transport properties (l.e.. Brownian diffusion coefficient, thermal diffusion 
factor) of small particles (considered simply as "heavy molecules") are 
calculated allowing for their variation with temperature across the boundary 
layer. The program dimensioning Is Increased to accurately obtain the sus- 
pended particle concentration profile Inside the much "thinner" (cf., thermal 
and momentum boundary layer thickness) mass transfer (Brownian diffusion) 
boundary layer for small particles (Sc » 1). 

STAN5 predictions of mass transfer coefficients are then used to examine 
the behavior and accuracy of correlations we recently developed to predict the 
enhancement In deposition rates due to thermophoresis In the presence of 
variable properties, transpiration cooling and/or viscous dissipation. These 
correlations, which successfully described results for self-similar laminar 
boundary layers, and law-of-the-wall turbulent boundary layers (ref. 6), are 
here found to be quite satisfactory and very economical, when applied locally 
to developing laminar and turbulent high Schmidt number mass transfer boundary 
layers . 

Whatever Is projected for electronic computational capabilities, more 
realistic system predictions/optimizations will Inevitably require subroutines 


7 See reference 18 for a comparison of STAN5 predictions with 
experimental heat transfer around a gas film-cooled cylinder In cross-flow. 
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which economically Incorporate the results of prior BL Integrations. In 
view of the complexity and the often prohibitive cost of running repetitive 
POE simulations for day-to-day engineering design calculations, the com- 
parisons discussed here suggest that local application of our previous mass 
transfer correlations (ref. 6) can cover typical cases of GT-Interest with 
acceptable accuracy. 
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FIG. 1 Accuracy of STAN 5 [11] and correlations [6, 7] in predicting 
the thermophoretic enhancementof particulate mass transfer rates 
when compared to self-similar laminar boundary layer calculations 
[11 for a flat plate; Sc m = 26, T w /Ta, = 0.6. 



FIG. 2 Comparison of correlations [6,7] with STAN 5 [111 predictions 
of mass transfer rate enhancement due to thermophoresis both with 
and without viscous dissipation; GT stator [15], Sc^ 26, T^T^ 0.6. 
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FIG. 3 Comparison of correlations [6.7] with STAN 5 [11] predictions 
of mass transfer rate enhancement due to thermophoresis; GT stator 
[15], Sc ro = 26, T w /T oo =0.& 



FIG. 4 Comparison of correlations [6,7] with STAN 5 [11] predictions 
of mass transfer rate enhancement due to thermophoresis; GT stator 
[15] , Sc^ = 100, T w /T ro - 0. a 
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FIG. 5 Comparison of correlations [6, 71 with STAN 5 
[111 predictions of mass transfer rate enhancement 
due to thermophoresis for a variable wall tempera- 
ture, GT stator [151 (yT* » 0.8— 0.5). 



FIG. 6 Comparison of correlations [6,71 with STAN 5 [111 predictions 
of mass transfer rate enhancement due to thermophoresis allowing 
transition to turbulence-, GT stator 1151, Sc OT = 26, T w IT ro = 0.5. 



FIG. 7 Comparison of correlations [6, 7J with STAN 5 (111 predictions 
of mass transfer rate reductions in the presence of gas transpiration 
cooling, both with and without thermophoresis; GT stator [15], 

Sc OT = 26, T 0 /Tod ” 0. 4. 
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